Kinetic Roughening in Slow Combustion of Paper 
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Results of experiments on the dynamics and kinetic rough- 
ening of one-dimensional slow-combustion fronts in three 
grades of paper are reported. Extensive averaging of the data 
allows a detailed analysis of the spatial and temporal devel- 
opment of the interface fluctuations. The asymptotic scaling 
properties, on long length and time scales, are well described 
by the Kardar-Parisi-Zhang (KPZ) equation with short-range, 
uncorrelated noise. To obtain a more detailed picture of the 
strong-coupling fixed point, characteristic of the KPZ uni- 
versality class, universal amplitude ratios, and the universal 
coupling constant are computed from the data and found to 
be in good agreement with theory. Below the spatial and tem- 
poral scales at which a cross-over takes place to the standard 
KPZ behavior, the fronts display higher apparent exponents 
and apparent multiscaling. In this regime the interface veloci- 
ties are spatially and temporally correlated, and the distribu- 
tion of the magnitudes of the effective noise has a power-law 
tail. The relation of the observed short-range behavior and 
the noise as determined from the local velocity fluctuations is 
discussed. 

PACS numbers: 64.60.Ht, 05.40.+j, 05.70.Ln 



I. INTRODUCTION 

The phenomenon of kinetic roughening of driven fronts 
is abundant in Nature. It has also become one of the 
paradigms in the physics of non-equilibrium systems. 
The theoretical interest arises since, as in critical equi- 
librium systems, the roughening behavior can be under- 
stood in terms of scaling exponents and scaling functions 
. A classification of kinetic roughening phenomena 
can be obtained by mapping the dynamics to an appro- 
priate Langevin equation which describes the interface 
dynamics, and also includes a noise term with system 
specific correlations and magnitude distribution. 

The simplest non-linear interface equation related to 
kinetic roughening is the celebrated Kardar-Parisi-Zhang 
(KPZ) equation 0, also related to the Burgers' equation 
of hydrodynamics. The KPZ equation for a growing sur- 
face can be expressed in the form 



dh(x,t) d 2 h{x,t) A 



dt 



dx z 



dh(x, t) 
dx 



+ F + V , (1) 



the case rj is short-range correlated in time and space, 
it can be substituted asymptotically by white, thermal 
(Gaussian) noise: (r](x,t)) = and (rj(x,t)r](x' ,t')) — 
2DS(x — x')S(t — t r ). The KPZ equation contains a linear 
diffusion term which imposes a surface tension through 
v, and the non-linearity proportional to A. The essen- 
tial physics of the KPZ universality class arises from the 
presence of additional symmetries, in particular a non- 
zero growth velocity component in the direction perpen- 
dicular to the local surface orientation. 

The scaling properties as implied by the KPZ equation 
can be discussed starting from e.g. the Family- Vicsek 
scaling H . It establishes a scaling relation for the surface 
width w 2 (L,t) = {(h~h) 2 ) 
system size L, such that 



as a function of time t and 



w 2 (t)=t 2f3 f(t/L z ) oc 



t 20 

L 2 * 



, for £ < L; 
, for £ > L. 



(2) 



Here £ ~ t x l z denotes a correlation length, which in- 
creases with time until the finite system size induces a 
saturation. The overbar and brackets denote spatial and 
noise averaging, respectively. Meanwhile three exponents 
have been defined: the growth (/3) and roughness (x) 
exponents, and the dynamical exponent z = x/P- F° r 
the KPZ equation with white noise ("thermal" KPZ or 
TKPZ), the exponents can be established exactly in one 
spatial dimension, and are (3 = 1/3, x — 1/2, and z = 3/2 



The TKPZ universality class is distinguished also by 
other characteristic properties, and not only by the scal- 
ing exponents. The scaling function f(x)of Eq. (|^) is 
another landmark of the KPZ equation In the 

steady state the KPZ interfaces are in fact equivalent to 
random walks, since a fluctuation-dissipation theorem ex- 
ists in (1 + 1) dimensions, depending only on the effective 
noise strength and surface tension, but not on A. On the 
other hand, temporal fluctuations on short time scales, 
and transient phenomena, do depend on the nonlinearity. 
Thus the scaling function f(x) and its properties provide 
a distinct confirmation of the TKPZ roughening. 

Consider the spatial [x 1 ) and temporal (t') two-point 
function of the interface fluctuations Sh(x, t) around the 
mean interface h(t) 



C 2 (r, t) = ([5h(x>, f ) - 8h(x> + r,t'+ t)] 2 ) a 



(3) 



where h(x, t) is the height of the interface, F the driving 
force, and r\ denotes the noise affecting the interface. In 



The equal-time limit t — t' defines G2M, an( l f° r t ne 
(T)KPZ universality class we have 



1 



G 2 (r) w Ar 2x . 



(4) 



Here A denotes an amplitude that depends on the 
parameters of the KPZ equation Eq. ([I]), A = D/u. Since 
this amplitude describes the steady-state fluctuations, it 
is related to the saturation width, and one finds §U^L^ 
that w sat = (A/12) 1 / 2 L X . Likewise, for the temporal 
correlation function, defined through Eq. (|J) by taking 
r = 0, we have IIhHlI 



(5) 



where = \\ 2 @ A° +1 Rq denotes a dynamical amplitude 
and includes a universal constant Rq that characterizes 
the universality class of the system. Notice that B de- 
pends on A, signalling that C s ^ is related to the transient 
interface properties. Rq can thus be determined as a uni- 
versal amplitude ratio, Rq = B / (\\\ 2 P A@ + ^). Other sim- 
ilar amplitude ratios can also be defined M. The values 
of Rq, as obtained |p|,|l0|-^2| from mode-coupling theory 
and simulations, are in the range Rg = 0.63 — 0.71. The 
correlation functions can as well be used to determine the 
universal coupling constant B such that 



9 =2 



.4 



B z/2 



1/X 



(6) 



Hwa and Frey have obtained || the value g* ~ 0.87 based 
on a mode-coupling solution of the (1 + l)-dimcnsional 
TKPZ problem. This calculation was later revised, and 
put in a more systematic form jl3| (see also ||), and 
compares well with the results from numerical studies of 
the scaling function ||l2| , [L0[ . 

Thus there are two distinct goals for an experimental 
study of interfaces, expected to be in the TKPZ univer- 
sality class. The power-law behaviors of the bare inter- 
face correlations should enable an estimation of the expo- 
nents x, (3 and z. Then, by augmenting these data with 
a determination of A, one can probe the scaling behavior 
beyond the scaling exponents by considering the univer- 
sal amplitude ratios and, possibly, the scaling functions 
of the fluctuations. Experimentally, the issue is compli- 
cated by the existence of cross-overs from short-range be- 
haviors to the asymptotic TKPZ behavior, which would 
affect e.g. the behavior of the height distribution, as 
compared to the pure TKPZ case. We have previously re- 
ported results on both the asymptotic scaling Jl^ l and on 
the properties of the short-range (SR) regime [|15| . In this 
paper our goal is to present a coherent picture of TKPZ 
physics in an experimental system, including more data 
and discussion related to the behavior at short scales. To 
this end, we have carried out additional high-resolution 
combustion-front experiments using very thin and light 
lens paper similar to that used in ]16[ , in addition to our 
previous data |l5| ] on two grades of ordinary copier pa- 
per, which are much thicker and heavier. The existence 
of large quantities of data allows us to do extensive av- 
eraging, typically over 10000 — 40000 interfaces. This is 



necessary to establish the correlation functions and their 
scaling properties reliably, in complete analogy with nu- 
merical experiments. Our results establish, clearly, the 
existence of asymptotic TKPZ scaling in terms of the 
exponents. Moreover, we are able to determine the uni- 
versal amplitude ratios and the fixed-point coupling con- 
stant in fair agreement with theoretical estimates. 

The short-range behavior in an experimental system 
can differ from the asymptotic scaling because of three 
fundamental reasons. First, the effective interface equa- 
tion may become TKPZ-like only during coarse-graining. 
A case of this kind is provided by the Kuramoto- 
Sivashinsky equation (T^j, where TKPZ scaling can be 
found only at large enough scales. Second, irrelevant 
terms in the renormalization group sense, neglected in 
the KPZ equation, may produce jl8| short-range effects. 
The third possibility is that, still within the experimen- 
tally accessible window of spatial and temporal scales, 
the effecive noise r\ in Eq. (|lj) proves to be more compli- 
cated. Consider a slowly moving interface, in the proxim- 
ity of a pinning/depinning transition that takes place at a 
critical F c . Since the interface moves in a landscape with 
frozen-in or quenched inhomogeneities, 77 will depend on 
h(x,t). For large enough driving forces, F » F c , such 
noise correlations vanish. For intermediate F one can de- 
fine a cross-over length scale, above which the behavior 
and noise correlations are TKPZ-like. Below that scale, 
very close to F c , the moving interface should have an 
effective (3 ~ \ — 0-75 p9| |, accompanied by apparent 
multiscaling. 

The correlations that arise from quenched noise can 
be generalized to temporal or spatial power-law correla- 
tions in the rj [f7pC|] with the correlator (r](x, t)rj{x' , t')) = 
2Dp(x~x')A(t—F). If the correlators p or A decay slowly 
enough, i.e. algebraically, the TKPZ exponents change 
0. The presence of cutoffs in the algebraic decay of p 
and A will induce again a cross-over to the TKPZ scaling. 
Similarly, a power-law amplitude probability distribution 
function P{rf) ~ > 1, can change the scaling 

behavior. For pure power-laws, there is a critical value 
of p c ~ 4, below which the roughening behavior depends 
on p. For p = 3, e.g., the q th order spatial two-point 
correlation functions have ||15| short range roughness ex- 
ponents that change from about 0.76 (q = 1) to 0.15 
(q = 9), while the asymptotic behavior is self-affine, with 
X- 0.76 p. 

In the following we will therefore also elucidate the SR 
behavior, by looking at the noise in the slow-combustion 
experiments, and how the correlations of this noise and 
its amplitude distribution are reflected in the SR scaling. 
Many different experiments on surface growth, erosion, 
imbibition (paper wetting), fluid invasion, and also slow- 
combustion fronts in paper, have failed to demonstrate 
asymptotic TKPZ scaling, though one can of course ask 
whether all of these problems should be described by a 
local interface equation p2j . Exceptional fluctuations in 
the noise amplitude have been observed in fluid-flow ex- 
periments p3[, and were reported earlier by us flq] for 
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slow-combustion fronts. The roles of the crossover, also 
noted in fluid-flow experiments |24f| , and in penetrating 
magnetic flux fronts in thin- film superconductors p5| , 
and of the effective noise in confusing the issue are to 
be noted. This may suffice as an explanation as to why 
TKPZ behavior was not seen in the first similar exper- 
iment on slow-combustion front roughening ]Ls| in con- 
trast with our results. 

Our impression is that in the slow-combustion fronts 
the SR physics largely arises from the presence of "non- 
ideal" noise related to the quenched impurities in the 
burning material, and to the natural fluctuations due e.g. 
to gas flow during the combustion process. We observe 
non-self- affine SR scaling, or multiscaling: the exponents 
X and (3 depend on the order of the correlation function 
from which they are measured. A high local front ve- 
locity is naturally coupled to a local, steep interface gra- 
dient. At such locations the interface propagates faster 
both parallel and perpendicular to the interface orien- 
tation, since A is positive. In the SR regime the front 
dynamics is coupled to the non-standard noise, which is 
both SR correlated and displays non-trivial amplitude 
fluctuations. The interplay of usual KPZ physics with 
the noise properties results in the absence of usual self- 
affine scaling. One of the origins of the quenched SR 
noise may lie in the local concentration variations of the 
potassium nitrate that is used to control the slow com- 
bustion. These experiments have been done at concen- 
trations large enough to ensure that an eventual pinning 
regime is far away (F ^> F c ). 

The paper begins in Section II with a detailed descrip- 
tion of the high-resolution experimental setup. Section 
III is devoted to an analysis of the paper samples on 
which the experiments are done. In Section IV we present 
results for the scaling exponents for all three paper grades 
used in the experiments, starting with the asymptotic 
long-range limit and touching then on the short-range 
behavior. Noise correlations are studied in detail for the 
lens paper, and compared with those for the copier papers 
reported in Ref. |]l5[ . Finally, we determine the univer- 
sal amplitude ratios for all three experimental cases in 
the asymptotic regime. Section V comprises a discussion 
about implications of our results to theory and to other 
related experiments. 

II. EXPERIMENTAL SETUP 

The experimental setup shown as a schematic diagram 
in Fig. |l] consists of a combustion chamber, a camera 
system, and a computer with necessary hardware and 
software. One side of the combustion chamber is made 
of glass and the opposite side is a detachable aluminium 
plate for installing the paper samples. The rest of the 
chamber sides are lined with a layer of porous material 
50 mm apart from the walls for making the incoming air 
flow laminar. In the middle of the chamber there is a de- 



tachable sample holder, designed for a maximum paper 
size of 600 x 400 mm 2 . It can be rotated with respect 
to the adjustable air flow so that convective transfer of 
heat ahead of the propagation front can be regulated, 
and the other flow-dependent features can be optimized. 
The sample holder is an open metallic frame whose sides 
are both lined with needles that keep the paper sheet 
planar during combustion. If the extra heat losses at the 
boundaries need to be compensated for, the sides of the 
sample can be heated with filaments that follow the com- 
bustion front. The volume flow through the chamber was 
adjusted if necessary with an electrical fan, placed on top 
of the chamber. The air flow through the chamber was 
also numerically simulated, and the velocity field around 
the sample was found to be quite homogeneous. To min- 
imize the dissipation of heat from the combustion front, 
and to keep the air flow laminar, we used relatively low 
volume flows (approximately 20 1/s for the copier papers 
and no forced flow for the lens paper). 

In the case the direction of air flow in the chamber was 
from bottom to top, combustion fronts were ignited from 
the top end of the paper samples by a tungsten heating 
wire in order to minimize convective heat transfer. In 
the original setup reported in Jli) , the recording of the 
front was done by a single charge-coupled device (CCD) 
camera, and the video signal was recorded on a Super 
VHS recorder. In the present [|l5| setup, the propagating 
front was recorded with three parallel PULNiX TM-6EX 
black and white CCD cameras, whose images are com- 
posed of 768 x 548 pixels. The largest combined image 
contains 2304 x 548 pixels. Cameras were attached to 
a stand that was operated with a pneumatic cylinder. 
The camera system was automatically moved in regular 
intervals along the direction of propagation. 

The digital frames of the three cameras were joined 
together and compressed online, and then recorded on a 
hard disk using a multi-level gray scale. The compressing 
of the frames was done by recording only a narrow stripe 
around the front line. Since recording was made in dark- 
ness, the only visible object was the combustion front. By 
omitting the dark backround, we were able to reduce the 
size of the individual frame file from 431 kB to approxi- 
mately 15 kB. The front height function was determined 
from each frame by first finding the pixels brighter than 
a given gray-scale value. A single-valued front line was 
fitted into the brightness profile in the strip found. In 
the intervals, where a front line could not be identified, a 
straight line was fitted by interpolation. The cylindrical 
image distortions caused by the lenses were corrected us- 
ing nonlinear warping. The method needs a collection of 
2D landmark points whose true locations are known to- 
gether with their distorted images. These were then used 
to define a global warping function. Using this warping 
function, corrections were made into the position data of 
the individual cameras before joining them. A series of 
typical fronts from the experiments are shown in Fig. (||). 

Improvements in data acquisition and introduction of 
the three-camera setup increased the temporal and spa- 
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tial resolution of the equipment. For a 390 mm wide and 
500 mm long paper sample, when the distance between 
the sample and the cameras was 70 mm, the instanta- 
neous recorded area was 310 x 74 mm 2 , and the pixel 
size was 0.135 mm. This was an order of magnitude be- 
low the maximum typical length scale related to areal 
mass variations in the paper samples, a few millimeters. 
It was also of the same order of magnitude, actually a bit 
smaller than the average width of the burning area in the 
combustion fronts, which sets the lower limit to details 
that can be detected. The maximum scanning frequency 
of the camera system was ten frames per second. Since 
the average velocity of the fronts for the copier paper used 
in the present experiments was approximately 0.5 mm/s, 
it was sufficient to store only every second frame of the 
digitized fronts. During one time step, At — 0.2 s, the 
front propagated less than one pixel length. For the lens 
paper we had to slightly modify our experimental setup. 
Needles in paper holder were not suitable for attaching 
the thin paper, instead metallic bars were used on both 
sides. Combustion fronts were ignited from the bottom 
and every successive front was recorded because of their 
much higher average velocity (6 — 8 mm/s) compared to 
those in copier papers, nd, as noted already above, for 
the lens-paper samples we used no forced ventilation in 
the chamber. This kind of setup was similar to that in 
the early experiments reported in |l6| . 

III. SAMPLES 

The experiments were done using two different, easily 
available grades of copier paper, with basis weights of 70 
gnr 2 and 80 gm~ 2 , and very thin lens paper with a basis 
weight of 9.1 gm~ 2 p6fl . It is important to realize that the 
composition of the paper has an impact on the propaga- 
tion of the combustion front. Slow-combustion fronts do 
not easily propagate in a material made of pure cellulose 
fibers only, because of their very low conductivity of heat. 
In normal copier papers approximately one fourth of the 
basis weight is composed of fillers. In our case, the fillers 
in both grades of copier paper were mostly calcium car- 
bonate (CaCOa) with a high heat capacity, in comparison 
with cellulose, which improved the propagation of com- 
bustion fronts. In addition, potassium nitrate (KNO3) 
was added as an oxygen source to all grades of paper to 
ensure a uniform propagation of slow-combustion fronts. 
The concentration of KNO3 in the samples was kept at a 
value of approximately 0.8 gm~ 2 . Potassium nitrate was 
usually added by embedding the samples in an aqueous 
solution of KNO3 for five minutes, after which they were 
dried in a press to maintain their planar shape. Due to 
the lower basis weight and much more porous structure 
of lens paper, drying of KNO3 solution was not homo- 
geneous over the sample, leading to an inhomogeneous 
concentration distribution of KNO3. This would have a 
significant effect on the experiments, as will be discussed 



later. To avoid such problems, KNO3 was added to lens 
paper samples by spraying. 

A second potential problem with paper samples is that 
it has been shown |27]] that, especially for low basis- 
weight laboratory paper sheets, there may exist nontriv- 
ial power-law correlations in the basis wight that may ex- 
tend up to about 15 times the fiber length, i.e. into the 
centimeter range. As it is well known |28j] that correlated 
noise affects the scaling behavior of kinetic roughening, 
maps of the local mass variations were prepared for both 
grades of copier paper before and after the KNO3 treat- 
ment. These (3 radiographs were taken from paper sam- 
ples of size 170 x 90 mm 2 using spatial resolution of 0.04 
mm. For the lens paper we used optically scanned im- 
ages to estimate the local mass variations. In addition, 
we measured the calcium and potassium concentration 
distributions in several layers of the same samples us- 
ing the laser- ablation method |^| with a spatial depth 
resolution of 0.2 mm. 

The radiographs of the copier paper samples were anal- 
ysed using the two-point density fluctuation correlation 
function 

C m (f) = {[m{x) — fn\[m(x + r) — fa]), (7) 

where m(x) is the local areal mass or the local basis 
weight, fh its average, and brackets denote an average 
over the sample. The results in Fig. |^ show that in both 
cases, after a distance of about one millimeter, a faster 
than algebraic (power-law) decay of correlations can be 
seen. This is in particular the case with the 70 gm -2 
paper, while in the heavier paper samples the noise level 
is reached rapidly. We may conclude that addition of 
KNO3 has not affected the basis weight correlations. The 
basis weight distributions P(m) = N{m) / N{m'), 
where N(m) is the number of locations with the same 
basis weight, differ from Gaussians only on the low basis- 
weight sides, which is most likely caused by rounding er- 
rors in the digitization of the (3 radiographs. 

The results of the laser-ablation measurements were 
analysed by a correlation function similar to that in 
Eq. (0). The difference was that, instead of the local 
basis weight, we used cither the calsium or potassium 
atomic emission line intensity, which is proportional to 
the vaporized mass. The intensity-fluctuation correla- 
tion functions were found to collapse to the noise level 
within a distance of one millimeter in every layer. Af- 
ter the KNO3 treatment the mean potassium intensity 
was found to be greater for the 70 gm~ 2 than for the 80 
g m -2 p a p er j n order to reach the same concentration 
of potassium nitrate for both grades, we used a stronger 
aqueous solution of KNO3 for the heavier paper. 

To summarize, the results indicate that there are no 
correlations in the filler, the local basis weight, or potas- 
sium nitrate distributions beyond a few millimeters, and 
that the structures of the samples are isotropic. 
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IV. KINETIC ROUGHENING OF SLOW 
COMBUSTION FRONTS 

A. Front width and correlation functions 

The most straightforward way to estimate the scal- 
ing exponents j3 and x f° r self-amne fronts is to use the 
scaling properties of the interface width w(L, t), as men- 
tioned in the Introduction. However, due to the large 
fluctuations in data when dealing with slow combustion 
fronts, it is useful to consider the q th order two-point 
height-difference correlation functions 



C q (r, t) = ([5h(x>, t>) - 5h(x> + r,t' + t)]") x >, t ', (8) 

where, again, h(x, t) is the height of the front at point 
x and time t, and 5h(x,t) = h(x,t) — h(t), and the bar 
denotes an average over a front while the brackets de- 
note an average over all configurations (fronts and burns) . 
Through this quantity, one can define the two functions 



and 



G g (r)=G„(r,0)~r**, 



(9) 



(10) 



which thus provide estimates for the roughness and 
growth exponents. In the saturated regime the functions 
G q (r) can be averaged over all times (steady-state con- 
figurations), and C Siq (t) over all spatial points. 

Another quantity of interest for which fluctuations can 
be efficiently averaged out, is the local width of the prop- 
agating front w(£,t), defined as 



w 2 (£,t) = (({h(x,t) - (h(x,t)) e ] 2 ) e ), 



(11) 



where the notation ()( denotes spatial averaging over all 
subsystems of size £ of a system of total size L. For 
growing self-affine interfaces, the scaling exponents sat- 
isfy (3 q = q(3 and \q = 9Xj an d l ne local width follows 
the Family- Vicsek scaling relation M given by 



w 2 (£,t) 



t 2li , for t < £ z 
£ 2 * , for t > £ z 



(12) 



This provides another way to estimate the scaling expo- 
nents from experimental data. 



B. Results for scaling exponents 

1. Data analysis 

In order to reduce the influence of boundary effects on 
the data, an area of width 270 mm taken in the mid- 
dle of the recorded area of width 390 mm was used to 
calculate h(x,t). The scaling exponents Xq an d (3 q were 
then determined by performing a linear least-squares fit 



to the corresponding two-point height correlation func- 
tions (Eq. (||)) in the scaling regime. An independent es- 
timate for the roughness exponent x was obtained from 
the local width. The early-time behavior of the surface 
width (Eq. ([ll])) also gave an estimate for the growth 
exponent (3. Accurate determination of the scaling ex- 
ponents was complicated by the "intrinsic" width of the 
fronts caused also by random structural inhomogeneities 
of the samples. Therefore, in the spirit of the usual con- 
volution ansatz []30"|-j32"[ by which a random process af- 
fecting the front but independent of the front dynam- 
ics will induce an additive constant in the square of the 
front width, we performed linear least-squares analyses 
of log 10 [G g (2r) - G q (r)] = qx^og w (r) + const, and of 
log 10 [C s ,g(2t) - C s , q (t)] = q/3\og w (t) + const, to get rid 
of additive constant factors in the correlation functions. 
The results given in the text are obtained after these "in- 
trinsic widths" were subtracted. 



2. Roughness exponent \ 

First, we present results for the roughness exponent xi 
using the second-order correlation functions and the local 
width. The main results for both y = X2/2 and (3 = /3 2 /2 
are summarized in Tables Q and y| for all three paper 
grades. In the latter Table all the fits have been done 
by first subtracting the intrisic widths from the data. 
The corresponding spatial correlation functions G 2 (r) are 
shown in Figs. ||and || so that in the latter figure the 
intrinsic widths have been subtracted. It is immediately 
evident from the data that there are two regimes of ap- 
parent scaling, separated by a crossover length r c . In the 
SR regime below the crossover length, G2 (r) scales with 
a rather large effective exponent of xsr — 0.90 for the 
heavier paper grades, and 0.85 for the lens paper. In the 
long-range (LR) regime r 3> r c , the effective exponent 
Xlr converges very close to the exact KPZ value of 1/2 
for all three cases. Beyond about 100 mm the statistics 
becomes worse and fluctuations larger. The inset in Fig. 
|B| shows estimates for the effective running exponents de- 
fined as Xcff = (1/2) log 10 [(G 2 (2r 2 ) - G 2 (r 2 ))/(G 2 (2n) - 
G 2 (ri))]/log 10 (r 2 /ri), where r 2 =4r*i. In all three cases, 
the LR asymptotic approach towards the TKPZ behav- 
ior is evident. It is also evident from this inset that no 
clear scaling regime can be found at short range with 
the higher apparent exponent. The behavior of the local 
width w 2 (l) averaged over the same data is consistent 
with G 2 (r) as is demonstrated in Fig. ^| where the intrin- 
sic width is subtracted. 



3. Growth exponent f5 

The temporal correlation functions G 2jS (t) for the three 
cases are shown in Fig. with the intrinsic widths re- 
moved. The poor time resolution (At =4.2 s) in the 70 
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gm -2 paper data prevented us from analysing in this case 
the temporal behavior very accurately. In the long-time 
regime the scaling behavior of the paper grades is again 
in agreement with the TKPZ case. The measured (3lr is, 
in particular after the intrinsic width has been removed, 
in good agreement with the TKPZ value 1/3. The largest 
deviations are shown by the lens paper as in this case the 
saturated regime was too short for an accurate determi- 
nation of (3. As was done in the case of the x's, we also 
used the local width, Eq. (|TT1), to compute the Plr- The 
scaling range is in this case slightly less than a decade, 
and the exponents obtained (without subtracting the in- 
trinsic width) agree well with the TKPZ one. Note that 
the same data for w(£, t) is used below to study its scaling 
function. 

In the data, there is a crossover at time t c from short- 
time behavior to asymptotic long-time regime visible in 
each case. The short-time growth exponent /?sr is again 
rather large, /3 SR = 0.75(5) for the 80 gm' 2 case, and 
0.64(3) for the lens paper. In the long-time regime be- 
yond the crossover time, the exponent /3lr is again fully 
consistent with the TKPZ value of 1/3. The inset in 
Fig. [?] shows the approach of the effective running expo- 
nents/3 cff = (1/2) log 10 [(C 2 , s (2i 2 )-C 2 , s (t 2 ))/(C , 2, s (2t 1 )- 
C 2 , s (ti ))]/ log 10 (i 2 /ii), where t 2 = 4ti, towards the 
TKPZ limit. Also in this case there is no clear scaling 
regime below the crossover time, although an apparent 
scaling exponent can be determined. 



4- Multiscahng at short range 

Spatial higher-order correlation functions, Eq. (|^), 
were determined to check the possible multiscahng prop- 
erties of the combustion fronts for the lens paper, to ac- 
company the data published earlier |l5| ] for the copier pa- 
pers. The behavior of the higher-order correlation func- 
tions indicates that, in the TKPZ regime, the interfaces 
are self-affine. On the other hand, in the SR regime the 
slopes of the spatial and temporal correlation functions 
of Figs. H depend on the order q of the correlation func- 
tion, i.e. show apparent multiscahng |33 34 2l]], similar 
to that found for the two grades of copier paper [ Ejl . 

Analysis of the fitted SR scaling exponents reveals that 
the xS-r(9)' s approach the value 1/2 for large q which 
means that larger local slopes follow a Gaussian distri- 
bution and are uncorrelated as in the long-range TKPZ 
regime. At the same time, PsrIiYs become very small 
and seem to approach zero. We conclude that the overall 
behavior of the SR exponents - with some uncertainty 
as for the higher momenta of the distributions - is very 
similar for all paper grades, and would thus seem to be 
a generic feature in these systems. One should however 
realize that uncorrelated random factors also appear in 
experiments, which affect the fronts and especially their 
scaling properties at SR if they involve high gradient val- 
ues. They are not related to front dynamics and appear 



therefore as "artefacts". Their effects can be analyzed 
qualitatively by using standard filtering techniques (e.g. 
median filtering). 



C. Noise correlations 

Next we consider the effective noise at the combustion 
fronts as determined from the fluctuations in the front 
velocity similarly to Ref. [^3). The physics behind this 
idea is that the presence of local large velocities can be 
related to the bare noise properties as the maximum front 
slopes are limited. The effective noise can be analyzed 
by considering the fluctuations of local velocities, defined 
as 5u(x, t) = u(x, t) — u(t), with u(x, t) = (l/r)[h(x, t + 
t) — h(x, t)], which obviously depend on the timescale r. 
Because of the crossover behavior, one has to analyze the 
noise correlations below (r < t c ) and above (r > t c ) the 
crossover scale. An alternative formulation for the noise 
amplitude is 



rj{x, t) = 8h(x, t + t) — Sh(x, t), 



(13) 



where 5h = h(x, t) — h(t). We have used this quantity to 
estimate the amplitude distribution, 



1 



(14) 



in the steady-state regime. Here N(r)) is the number 
of positions where r\ has the same value, and the sum 
goes over all values of r\. We determined rj(x,t) for time 
steps r = 0.5,1.0,2.0,4.0, and 8.0 s for the lens paper, 
and show the related distributions P(rj) for -q > in 
Fig. [I For short time steps the distributions have to 
a good degree of accuracy a power-law tail of the form 
P{rj) = crj~(>- l+1 \ with fi ~ 1.7. For increasing r the 
power-law contribution in the tail of P(j]) becomes less 
visible, the exponent /i increases towards /i = 5, and 
the distribution approaches a Gaussian. The behavior is 
similar to that found for the copier paper ]l5| , except for 
the exact value of the exponent \x which was found to be 
about 2.7 at short times (r = 1.0 s). 

The local velocity fluctuations were calculated using 
time steps r = 0.5, 1.0, 2.0, and 4.0 s for the lens paper. 
In the inset of Fig. || we show the corresponding distri- 
bution function P(5u). For large t's, P(Su) seems to be 
nearly Gaussian, as expected. For small t's, it has an 
asymmetric tail towards higher velocities, which is con- 
sistent with nonzero skewness of the TKPZ height distri- 
bution [Q. We note that ash formation increases local 
heat transfer ahead of the combustion front, contributing 
to this effect in the velocity distribution. 

The two-point correlation function of these fluctua- 
tions can be expressed in the form 



C u (x, t) = {5u(x + x, t + t)5u(x , t )). 



(15) 



G 



This expression is useful for determining possible spatial 
and temporal correlations. In Fig. [lo] we show the spa- 
tial correlation functions C u (x, 0) for different time steps 
t, and the corresponding temporal correlation functions 
C u (0,t) in the inset. The effective spatial and temporal 
noise was found to be uncorrelated above the crossover 
scales r c and t c , which is in agreement with the TKPZ 
behavior. Also in this case the behavior of lens paper 
is analogous to that of the copier papers for which data 
have been reported in Ref . |ll| . Notice that the crossover 
time is longer than the ratio of the spatial crossover scale 
and the front velocity would predict. 



D. Universal amplitude ratios 

We now look at such properties of the TKPZ class as 
the universal amplitude ratios and the universal coupling 
constant, discussed in the Introduction. The experimen- 
tal determination of these quantities is slightly hampered 
again by the presence of the crossover scales to the TKPZ 
behavior. Thus the correlation functions and the nonlin- 
earity A need to be considered in the asymptotic regime 
only. This means that the available scaling range is lim- 
ited, but, as we shall see, the results are still in rather 
good agreement with the theoretical predictions. In any 
case, correlation functions provide more accurate values 
for these quantities than the front widths, and we there- 
fore use them here. 

We measure from the spatial and temporal correlation 
functions the prefactors A and B (see the Introduction 
for a discussion). Figure |ll| demonstrates the procedure 
by which A and B were determined in the saturated 
TKPZ regime for all three paper grades using reason- 
able choices for the saturated parts of the correlation 
functions. The actual values for A, B, and the other pa- 
rameters are listed in Table III. The amplitude values 



were augmented by an independent measurement of A. 
The interface velocity was determined at each location 
with a time interval r long enough to attain the TKPZ 
regime, with a simultaneous determination of the local 
interface slope. The velocity vs. slope data were then 
fitted by a parabola using slopes between -0.5 and 0.5 
p5[ , corresponding typically to about a half of all the 
slopes. 

Combining the results for A, A, and B with the exact 
TKPZ exponents, we obtain for the coupling constant g* 
the values shown in Table III . It is evident from this table 



that the results are in good agreement with theory, with 
the largest deviation showing up for the 70 gm~ 2 grade. 
Recall that the mode-coupling value is g* ~ 0.87 [||. 
If we use the measured exponents instead, g* becomes 
0.8(2), 0.8(2), and 2.7(4) for the 70 gm~ 2 , 80 gm~ 2 and 
9.1 gm~ 2 grades, respectively. The last value is due to 
the large effective (3 for the lens paper, caused by the 
very short saturated regime, and thereby poor statistics 
in the temporal correlation function. The same quanti- 



ties, A, A, B, and the TKPZ exponents, also define the 
universal amplitude ratio Rq , whose measured values are 
also shown in Table III . As the mode-coupling and simu- 
lation results for Rq vary in the range Rg = 0.63 — 0.71, 
agreement with theoretical predictions is good also in 
this case. Notice that g* and Rg are functionally re- 
lated when the exact TKPZ exponents are used in their 
expressions. 

The local width w(£, t) allows furthermore for the 
determination of a scaling function related to that of 
Eq. (g). One can show, see e.g. that w(£,t) = 

Cp£ x F(\X\Cpt/ '£ z ), where Cp is determined by the asymp- 
totic value (at t = oo) of w, w(£) — Cp£ x . Using a maxi- 
mum value of £ ~ 14 cm, and fitting the Cp values for the 
three grades, one can collapse the time-dependent data 
with the aid of the A values as measured from the slope- 
dependent local velocities of the fronts |^C|pT| . The result 
is shown in Fig. |l2|. As the measured A and Cp values are 
not very accurate, we show in the inset of this figure the 
best collapse of the data achieved by using Cp essentially 
as a free parameter. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have made a thorough survey and 
analysis of our experiments on kinetic roughening of slow- 
combustion fronts in paper. The general conclusion is 
that well-controlled experiments lead to a clear asymp- 
totic scaling that is unequivocally in accord with the ther- 
mal KPZ universality class. This is evident from the 
asymptotic behavior of the gth order correlation func- 
tions, and from the scaling exponents determined from 
them, as well as from independent measurements us- 
ing the local width. We also augment the evidence by 
computing a universal amplitude ratio characterizing the 
TKPZ universality class, and the universal (KPZ fixed- 
point) coupling constant, and find good agreement with 
theoretical estimates. 

In the analysis of our experimental results it proved 
to be of crucial importance to do extensive averaging 
over independent slow-combustion fronts. Fluctuations 
in the noise affecting the fronts give rise to wide fluc- 
tuations in individual fronts, also from burn to burn. 
There is quenched noise due to density variations typical 
of paper-like materials, and the effective noise also in- 
cludes dynamical effects. A burning front creates around 
it a fluctuating flow of air, and this flow will give rise to 
fluctuations in the "effective" heat conductivity (due to 
convective transfer of heat) and in cooling. These dy- 
namical fluctuations could be demonstrated by adjust- 
ing the overall air flow in the combustion chamber, and 
thereby regulating the crossover scales parallel to the 
flow. The temporal crossover scales are always clearly 
longer than the spatial crossover scales along the fronts 
divided by the average front velocities. This is indicative 
of dynamical convective effects. We found that averaging 
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over approximately ten independent burns is needed for 
fairly reliable estimates for the quantities measured. For 
copier paper this means in practice averaging over about 
40000 individual fronts, and for lens paper, with much 
faster propagating fronts, over several thousand individ- 
ual fronts. For the lens paper, averaging was hampered 
by the shortness of the saturated regime. 

In addition to the demonstration of the asymptotic 
TKPZ behavior, another interesting feature is the persis- 
tent crossover present in slow combustion of paper. The 
apparent scaling properties of the fronts are markedly 
different at short time and length scales as compared to 
the asymptotics. This phenomenon was first observed 
[24 1 in fluid-flow experiments, and very recently also in 
[25| penetrating flux fronts in thin-film superconductors. 
As discussed in the Introduction, several possible mecha- 
nisms have been suggested for this phenomenon (see also 
Ref. 0] for a review). Before our work, the reasons and 
consequences of the two regimes was not settled experi- 
mentally though power-law distributed amplitudes in the 
effective noise were reported both in [^| for fluid-flow ex- 
periments and in [ plj| for slow combustion of paper. In 
the latter experiments, performed by us, noise correla- 
tions were furthermore found to be short ranged, both in 
space and in time. Moreover, the "decay lengths" of these 
correlations seem roughly to coincide with the crossover 
scales in the height-height correlation functions. It thus 
appears that the SR correlations in the noise affecting the 
fronts, either quenched or annealed, or both, are likely to 
be responsible for the crossovers and the related higher 
apparent exponents at short range. The continuous de- 
cay of SR correlations would also explain the lack of true 
scaling in that regime, i.e. the running exponents do 
not show plateaux there. In the case of slow combustion 
of paper at least, the effective noise clearly is partly of 
dynamical origin. Thus quenched noise alone cannot be 
decisive for the short-range effects, although it leads to a 
similar crossover behavior. Our measured SR exponents 
shown in Table [n] do not agree with the apparent SR 
exponents for moving fronts with quenched noise: close 
to the depinning transition /3 ~ ^ — 0.75 jfj|. Some ef- 
fects at SR due to the irrelevant terms in the appropriate 
Langevin equation cannot of course be ruled out. 

Multiscaling at short range seems, on the other hand, 
to result from other features in the fronts. It is evi- 
dent that regions of high gradient (in absolute value) are 
amplified in the higher-order correlation functions, and 
thereby affect the multiscaling properties of the fronts. 
Distinct regions of such high gradient values appear if 
successive fronts are plotted as surface diagrams. Note 
that digitizing errors, and sharp natural defects in the 
observed data, also make a contribution to the appar- 
ent multiscaling properties, and should be filtered out 
as artefacts. Fairly sharp steps, which occasionally ap- 
pear in the (filtered) fronts, begin to move also sideways, 
driven by the nonlinear term in the KPZ equation Eq. (Fy). 
These moving structures appear as local "avalanches" in 
the propagating fronts. By analyzing distinct intervals 



of propagating combustions fronts, containing varying 
amounts of such avalanches, we can conclude that the 
q-dependence at SR of the slopes of the spatial and tem- 
poral correlation functions extend longer when there are 
more avalanches present. This effect is however difficult 
to quantify precisely as it is not possible to completely 
remove the avalanches, nor the experimental artefacts, 
from real data. 

We have only considered here well-propagating fronts 
in which pinning effects do not appear. By reducing the 
KNO3 concentration, we can however approach the pin- 
ning regime in a rather controlled way. A lot of statistics 
must be produced again, since the exact amount of KNO3 
absorbed in the paper samples (the "driving force" ) is dif- 
ficult to regulate. Nevertheless, in this way we think one 
can eventually determine the interface behavior also at 
pinning. It is evident that there will also be crossovers 
present in this case, with the actual crossover scales prob- 
ably depending on the driving force. 
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FIG. 2. Series of typical digitised fronts. The time step 
between successive fronts is 10 s, and the width of the digitised 
area is 310 mm. 
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FIG. 3. A log-log plot of the correlation function C m (r) 
for four copier paper sheets of basis weights 70 gm~ 2 , 70 
gm -2 treated with KNO3, 80 gm~ 2 , and 80 gm -2 treated 
with KNO3 (from top to bottom). The dashed line shows 
Gaussian behavior. The curves have been shifted vertically 
for clarity. 



FIG. 1. A schematic diagram of the experimental setup. 
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FIG. 4. The spatial correlation function Gi(r) vs. r aver- 
aged over 10 burns of the 70 gm -2 (o) and 6 burns of the 80 
gm~ 2 copier paper (□), and 32 burns of the 9.1 gm -2 lens 
paper (★). The crossover lengths r c are 5.2(1), 7.9(2), and 
11.1(1) mm for these paper grades, respectively, shown with 
vertical lines, and the solid lines indicate slopes that corre- 
spond to the exponents in Table |. The curves have been 
shifted for clarity. 
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FIG. 6. The local width after the intrinsic width was 
subtracted [w 2 (2l) — iti 2 ^)] 1 / 2 vs. I averaged over 11 burns 



of the 70 gm 



and 18 burns of the 80 gm copier paper 



(□), and 24 burns of the 9.1 gm -2 lens paper (*). Solid lines 
indicate slopes that correspond to Table III The curves have 
been shifted for clarity. Inset shows the effective rougheness 
exponent \cfi- 
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FIG. 5. The spatial correlation function after the intrinsic 
width has been subtracted G2 (2r) — Gi (r) vs. r averaged over 
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FIG. 7. The temporal correlation function after the intrin- 
sic width was subtracted C s ,2(2t) — C a ,2(t) vs. t averaged over 



10 burns of the 70 gm (o) and 6 burns of the 80 gm copier 11 burns of the 70 gm (o) and 18 burns of the 80 gm 



paper (□), and 32 burns of the 9.1 gm - lens paper (*). The 
solid lines indicate slopes that correspond to the experiments 
in Table |[lj The curves have been shifted for clarity. Inset 
shows the effective rougheness exponent XcS vs. r, where the 
horizontal line is \ = 1/2. 



copier paper (□), and 24 burns of the 9.1 gm -2 lens paper 
(*). The curves have been shifted for clarity. The solid lines 
indicate slopes that correspond to the experiments in Table 
[n]. Inset shows the effective growth exponent /3 e ff vs. t. The 
horizontal line denotes /3 = 1/3. 
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FIG. 8. (a) The q th order spatial correlation function 



G q {r) 



1/9 



r averaged over 32 burns of the 9.1 gm 



lens paper, (b) The q th order temporal correlation function 
C a , q (t) 1 / q vs. t averaged over 24 burns of the 9.1 gm~ 2 lens 
paper. The functions in (a) and (b) are plotted for even mo- 
ments from q — 2 to 10. 
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FIG. 9. The noise amplitude distribution P(ri > 0) aver- 
aged over 35 burns of the 9.1 gm -2 lens paper with an av- 
erage velocity > 8.4 mm/s. The distribution was calculated 
for time intervals r = 0.5, 1.0, 2.0, 4.0, and 8.0 s shown in the 
figure. The linear fits have slopes -2.71, -2.67, -2.72, -2.99, 
and —4.97, left to right. The inset shows the local velocity 
fluctuation distributions for the same time intervals r. 
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FIG. 10. The spatial and temporal velocity fluctuation cor- 
relations C u {x,0) and C„(0,t) (inset), respectively, for time 
intervals r = 0.5, 1.0, 2.0, and 4.0 s for the 9.1 gm -2 lens pa- 
per. The crossover length r c and time t c (inset) are shown by 
vertical lines. 
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FIG. 11. (a) Measured amplitude of the spatial correlation 
function A = r~ 2x G2(r) vs. r averaged over 10 burns of 
the 70 gm -2 (o) and 5 burns of the 80 gm~ 2 copier paper 
(□), and 32 burns of the 9.1 gm -2 lens paper (*) (inset). 
The horizontal lines indicate the value estimated for A for 
different paper grades, and it has values 0.54, 0.42, and 3.47, 
respectively. The roughness exponent was fixed at \ = 1/2- 
(b) Measured amplitude of the temporal correlation function 
B = t~ 213 C s ,2(t) vs. t averaged over 11 burns of the 70 gm -2 
(o) and 18 burns of the 80 gm~ 2 copier paper (□), and 24 
burns of the 9.1 gm -2 lens paper (*) (inset). The horizontal 
lines indicate the value estimated for B for different paper 
grades, and it has values 0.18, 0.15, and 9.01, respectively. 
The growth exponent was fixed at /? = 1/3. 
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FIG. 12. Scaling plots for w(i,t) = £ 1/2 C e F (\\\C e t/t i/2 ) 
averaged over 39 burns of the 70 gm~ 2 (o) and 34 burns of the 
80 gm~ 2 copier paper (□), and 24 burns of the 9.1 gm~ 2 lens 
paper (*). The inset shows the best collapse of these data, 
achieved by letting Ce vary (a variation within a few percent 
was enough) as a free parameter. 
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TABLE I. The scaling exponents (5 and x, crossover values r c and t c , and average velocities. 





70 gm" 2 


C 2 (r,t) 
80 gm" 2 


9.1 gm" 2 


70 gm" 2 


w(£,t) 
80 gm -2 


9.1 gm" 2 


XSR 
XLR 


0.88(2) 
0.53(3) 


0.89(2) 
0.51(3) 


0.83(1) 
0.53(4) 


0.81(6) 
0.57(1) 


0.83(5) 
0.55(2) 


0.81(1) 
0.56(2) 


r c [mm] 


4.7(4) 


6.0(5) 


11(2) 


12(7) 


14(6) 


18.8(4) 


/3sr 

0LR 


0.59(4) 
0.40(3) 


0.69(2) 
0.39(3) 


0.61(2) 
0.46(2) 


0.29(3) 


0.32(3) 


0.28(5) 


tc [s] 


25(9) 


27(5) 


3.7(4) 








v [mm/s] 


0.64(2) 


0.50(2) 


9.0(5) 





TABLE II. The scaling exponents (3 and x obtained by first subtracting the intrinsic widths from the data. 





70 gm" 2 


C 2 (r,t) 
80 gm -2 


9.1 gm" 2 


70 gnr 2 


w(£,t) 
80 gm -2 


9.1 gm" 2 


XSR 

Xlr 


0.90(3) 
0.50(4) 


0.90(4) 
0.47(4) 


0.85(1) 
0.50(6) 


0.84(6) 
0.56(1) 


0.87(8) 
0.52(5) 


0.95(6) 
0.51(1) 


PSR 

Plr 


0.36(3) 


0.75(5) 
0.34(4) 


0.64(3) 
0.43(6) 









TABLE III. Results for the correlation amplitudes, A, and the universal quantities using f3 = 1/3 and \ — 1/2- 



Paper grade 


A 


B 


A 


Rg 


9 


70 gm" 2 


0.52(2) 


0.186(12) 


0.465(2) 


0.74(6) 


0.79(9) 


80 gm" 2 


0.475(7) 


0.14(1) 


0.370(1) 


0.73(5) 


0.76(8) 


9.1 gm" 2 


3.4(1) 


8.0(8) 


4.0(4) 


0.62(8) 


1.0(2) 
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